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\zi^- Abstract. We discuss some basic aspects of the dynamics of a homogenous 

^^^ • Fermi gas in a weak random potential, under negUgence of the particle pair in- 

teractions. We derive the kinetic scaling limit for the momentum distribution 
function with a translation invariant initial state and prove that it is deter- 
mined by a linear Boltzmann equation. Moreover, we prove that if the initial 
state is quasifree, then the time evolved state, averaged over the randomness, 
(^ ' has a quasifree kinetic limit. We show that the momentum distributions de- 

Q,, termined by the Gibbs states of a free fermion field are stationary solutions of 

I , the linear Boltzmann equation; this includes the limit of zero temperature. 



^ 1. Introduction 

m 

>. 

0^ ^ We investigate the Boltzmann limit for the dynamics of a quantized field of non- 

^^ • relativistic electrons in a disordered medium. Our approach is based on techniques 

^ ; developed in [TUl m [H [H] and [3 [5] , used to derive the Boltzmann limit for 

the one-particle Anderson model at weak disorders; see also [1TJ[TH]. We refer also 

to [ll[5l[i[18l[20] for related works. 



o 



jCr* ■ We consider a gas of fermions on the lattice A^ := [—■§■, ■|]'^ n Z*^ in dimension 

^—^ ' d > 3 and with periodic boundary conditions, for L ;» 1. We denote the dual 

lattice by A^ — ^l/L, and write J dp = jd X^neA* ^'^^ brevity. Letting ^ = 
^ , ®n>o Ai ^^(^i) denote the Fock space accounting for scalar fermions on A^, we 

JH I denote the creation- and annihilation operators by aj" , Op, with p e A^, satisfying 



C^ 



the usual canonical anticommutation relations. 

Let 21 denote the C*-algebra of bounded operators on ^. We let po denote a 
translation invariant, normalized state on 2t, which preserves the particle num- 
ber (i.e., po{NA) = po{AN) for all A € Qi, where N = J^x^'t'^x is the number 
operator). 

We consider the Hamiltonian 

H^ ■■= / dpE{p)a+ap + i]V^ (1.1) 

which generates the dynamics of a free Fermi gas coupled to a random potential 

K; := ^ Lo^a'^ax (1.2) 

aJSAi 

1 
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where {uJx}xgAl ^-re real Gaussian i.i.d. random variables, and < 77 ^ 1 is a 
small coupling constant accounting for the disorder strength. We assume that the 
kinetic energy function is given by 

d 

^(p) = ^cos(2^p,), (1.3) 

i=i 
i.e., the Fourier multiplication operator determined by the centered nearest neighbor 
Laplacian (A/)(x) = J2\y~x\=if(.y) on Z'^. 

We are interested in the long-time dynamics of the fcrmion field described by 

PtiA) -.^ poie^'"-Ae-^'"-), (1.4) 

where A e 21. While we are neglecting the pair interactions between the electrons, 
the effective interaction between the particles through their coupling to the random 
potential, and due to the Pauli principle remain significant. We prove the following. 
In a time scale t = ^ where T > denotes a macroscopic time variable, we find, 
in the thermodynamic limit, that for all T > and for all test functions /, g of 
Schwartz class iS(T'*), 

4')(/;g) := lim lim E[pt/^2( a+(/) 0(5) ) ] - / dp Ft{p)W) 9{p) . (1-5) 

where Ft{p) satisfies the linear Boltzmann equation 

drFrip) = 2 7r f du5{E{u) ^ E{p)) {Ft{u) - Ft{p)) (1.6) 



with initial condition Fo{p) = liniL^oo t7Po( OpOp)- The proof is based on a 
generalization of methods due to Erdos and Yau in 14 , and extended in [7], for the 
derivation of linear Boltzmann equations from the random Schrodinger dynamics 
in the weakly disordered 1-particle Anderson model. 

We observe that if po is the Gibbs distribution of the free fermion field, the 
corresponding momentum occupation density (the Fermi-Dirac distribution) 

for inverse temperature (3 and chemical potential /i, is a stationary solution of the 
linear Boltzmann equation (|1.6p . for all /3 > 0. This is also valid in the zero 
temperature limit /3 — > cxd where in the weak sense, 

1 + Je(,,)-,) ^x[i^(p)<M], (1.8) 

which is nontrivial if /x > 0. Erdos, Salmhofer and Yau have proved in their land- 
mark work [inilllllll] that for a time t beyond the kinetic scale 77"^, the effective 
dynamics of a single electron is diffusive; i.e., in this time scale, a wave packet 
evolves in position space according to the solution of a heat equation. Accord- 
ingly, we expect the Fermi-Dirac distribution to remain a stationary solution in the 
diffusive limit, and the corresponding time scale addressed in [lO l ITT | IT2j . 



We remark that the translation invariant model without the random potential 
(i.e., 77 = 0) but including the full repulsive particle pair interaction is determined 
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by the Hamiltonian 

Hx := j dp E{p) a'^ ap + \ ^ Uy a^ v{x - y) Uxtty . (1.9) 

It is widely believed that in a time scale t = -^, the momentum density Ft{p) '■— 
limA^o liniL^oo 77-PT/A2(ap Op ) for the dynamics generated by Hx satisfies the 
Boltzmann- Uhlenbeck- Uehling equation 

drFTip) = -An dpidp2dqidq2\v{pi- qi) -v{pi- q2)\'^ 5{p--pi) 
S{pi +P2-qi-q2) 5{E{pi) + E{p2) - E(qi) - E{q2) ) 
FT{pi)FT[p2)FT[qi)FT{q2) ~ FT[qi)FT{q2)FT{pi)FT{p2) 

where Ft{p) := 1 — Ft{p). The derivation of (jl.lOp from the microscopic quantum 
dynamics is an extremely challenging open problem; for some work in this direction, 
see [4j [131 HB [22] . We note that (|1.7p is also an equilibrium solution of (|1.10|) , which 
is ealily seen by noting that Fo(p) = e^'^^'^^'^" ^^ Ff){p) . This is a consequence of the 
fact that (|1.7p is a function of the kinetic energy E{p) which is a collision invariant 
in both (|1.6p and (|1.10p . As a matter of fact, any distribution of the form f{E{p)) 
is stationary for (|1.6p : on the other hand, the special structure of (|1.7p is necessary 
for it to be a stationary solution of (|1.10p . For a combined Boltzmann limit of the 
coupled model with A, 77 > (which is an open problem) we believe that the kinetic 
energy E{p) will remain a collision invariant, and that the momentum distribution 
(|1.7p will remain a stationary solution of the resulting Boltzmann equation. 

A contextually related question is the problem of the stability of the Fermi sea 
for a gas of interacting fermions. This is an important problem in mathematical 
physics which has in recent years received much attention, especially due to the 
landmark works of Feldman, Knorrer, and Trubowitz summarized in [15] . 

We also address the question how strongly the electrons are effectively correlated 
through their interactions with the random potential. To this end, we assume that 
Po is number preserving, homogenous, and quasifree. That is, for any tuple of test 
functions fi, . . . , fr, 9i, ■ ■ ■ ,gs, 

Paia+ifi)---a+{fr)a{gi)---a{gs)) = (5^,^ dct[po(a+(/j) a(gf) ) ]^^^^i . (1.10) 

We consider the dynamics generated by H^, and observe that since H^ is bilinear 
in a'^jU, the time evolved state pt is almost surely quasifree. However, the state 
averaged over the randomness is not quasifree, 

lim E[ptifi,...Jr;gi,...,gr)] ^ det[ lim E[pt(/,- ; ^f)] ]J^^i , (1.11) 

L — >oc L — 'oc 

for any i > if 77 > 0. This is not surprising because quasifreeness is a nonlinear 
condition. We prove that in the kinetic scaling limit stated above, the limiting 
2r-correlation functions are quasifree, 

^xih^-'-Jr; gi,...,gr) 

:= lini lim E[/?t/,,2 ( a+(/i) • ■ • a+(/r) a(gi) • • • a(.g2) ) ] 

T)^0 L—foo 

= det[d^\f,;g,)]l,^, (1.12) 
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for any r e N. The proof is based on an extension of the proof in [^ for the 1- 
particle Anderson model at weak disorders that the random Schrodinger evolution 
converges in arbitrary higher mean to a linear Boltzniann evolution. Quasifreeness 
of the 2r-point correlation functions is a significant ingredient in some approaches 
to the problem of quantum charge transport; see for instance [2 and the references 
therein. 



2. Definition of the model 

We give a detailed definition of the mathematical model described in the previous 
section. We consider a fermion gas in a finite box A^ := [— -j, ■|]''nZ'^ of side length 
L 3> 1, with periodic boundary conditions, in dimensions d > 3. We denote its dual 
lattice by A^ := A^/L C T''. For the Fourier transform, we use the convention 



where p G A^, and 



m := J2 ^"'"'" /(^) ' (2-1) 

xEAl 

/(^) = i^ E e'^^'-^/b) (2-2) 

for its inverse. For brevity, we will use the notation 

fdp^^ii: (2-3) 

in the sequel, which recovers its usual meaning in the thermodynamic limit L -^ cxd. 
We denote the fermionic Fock space of scalar electrons by 

5 = ;?« , (2.4) 

n>0 

where 

n 

So = C , dn = f\i\^L), n>\. (2.5) 

1 

We introduce creation- and annihilation operators a^ , a^, for p, q (z A^, satisfying 
the canonical anticommutation relations 

a+ a, + a, a+ = <5(p - -z) :- { q fth'Irwis'e. (^.6) 

We define the fermionic manybody Hamiltonian 

H^ ■.= T + r,V^ (2.7) 

where 

T = f dpE{p)a+ap (2.8) 

is the kinetic energy operator, and 

VL := ^ ujxat'^x (2.9) 

x£Al 
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couples the fcrmions to a static random potential; {uJx}xeAL i^ ^ field of i.i.d. real- 
valued random variables which wc assume to be centered, normalized, and Gaussian 
for simplicity. Thus, 

E[c^,] = , E[luI] = 1 (2.10) 

for all X G A^^. Moreover, we assume that 

d 

E{p) ^ Y.^os{2t:pj) , (2.11) 

which defines the Fourier multiplier corresponding to the nearest neighbor Laplacian 

onZ'^. 

Let 

N -.^ Y. «-°- (2.12) 

denote the particle number operator. It is clear that 

[H^,N]^0 (2.13) 

holds. 

Let 21 denote the C* -algebra of bounded operators on 'S- We consider the dy- 
namics on 21 given by 

at{A) ^ e'*"-^ Ae-'*"- (2.14) 

generated by the random Hamiltonian H^. 

3. Statement of the main results 

We consider a normalized, translation- invariant, deterministic state 

po:2l — >C. (3.1) 

We define the time-evolved state 

Pt{A) := po{e^'"-Aer^'"-), (3.2) 

with t G R, and initial condition given by pQ. We particularly focus on the dynamics 
of the averaged two-point functions 

nptia^a,)], (3.3) 

where p,q E A^ . Clearly, 

E[po{a+aq)] ^ po{ajaq) ^ 5{p - q) j^ pn{a+a.p) , (3.4) 

where 

8{k) :- L'^4, (3.5) 

and where 



1 li p ^ q 
otherwise 



(3.6) 
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denotes the Kronecker delta on the lattice A^ (mod T ). We remark that for 
fermions, 



1 

L' 
since \\ap \\ — L'''^ in operator norm, Vp G A^ 



< ;^Po(<ap) < 1, (3.7) 



3.1. The Boltzmann limit. We denote the microscopic time, position, and veloc- 
ity variables by (t, x, p) , and the corresponding macroscopic variables by (T, X, V) = 
{r]'^t,r]'^x,v). We prove that the momentum distribution ft{q) converges to a solu- 
tion of a linear Boltzmann equation in the limit 77 ^ 0. 

Theorem 3.1. We assume that po is translation invariant. Then, the averaged 
two-point functions are translation invariant, 



E[pt{a+if)a{g))] = J dp f{p) g{p)nptia;ap)] , (3.8) 

(i.e., diagonal in a+, Op) for any /, g G (S(T'') of Schwartz class, and the thermody- 
namic limit 

^p''\f:9) ■■= lim nPT/n4a+{f)a{g))] (3.9) 

exists for all f,g^ S{T'^), and T > 0. 

For any T > and all f,g G S{T'^), the limit 

nP{f-g) := Yim^p''\f;g) (3.10) 

exists, and is the inner product of f,g with respect to a Borel measure FT{p)dp, 

^P(f-,g) = j dpFT{p)W)g{p), (3.11) 

where Ft{V) satisfies the linear Boltzmann equation 

drpFriV) = 2n f dU 5{E{U) - E{V)){Ft{U) - Ft{V)) , (3.12) 

with initial condition 

Fo{p) = lim 7^Po(ap.. ap^* ) (3-13) 

forp&T^, where PA* := Qj, (p) fl A^, and Qs{p) -^ p + [—S,6y. 

We note that there exists a unique pA* G A|^ such that \p — pA* \ < ^, for every 

peT"^. 

An initial condition of particular interest is the Gibbs state (with inverse tem- 
perature (3 and chemical potential p) for a non-interacting fcrmion gas, 

Po{A) = ^Tr(e-'^(^-^^)A) (3.14) 

where ^/3,/j ■— Tr( e^^*-"^^^^' ). The corresponding momentum distribution func- 
tion 

£m _p„(«+ap) = ^^^^(^(^)_^) (3.15) 
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is a fixed point of the linear Boltzmann equation (I3.12|) . for all /? > 0, including the 
zero temperature limit (3 -^ oo where in the weak sense, 

1 + Je,,,-,, - X[E{P) < M] , (3.16) 

which is nontrivial if /i > 0. We note that all our results in this paper remain valid 
in the limit f3 —>■ oo. 

3.2. Quasifreeness. We prove that if in addition to the conditions formulated 
above, the initial state po is quasifree, then E[pt], which is not quasifree for any 
t > if 77 > 0, becomes quasifree in the kinetic scaling limit of Theorem l3.1l 

A state pq is quasifree if for any normal ordered product of creation- and anni- 
hilation operators 

«Pi ••■"p,-"«i ■■•«9= ' (3-17) 

with arbitrary r, s G N and pi, Qj G A^, 

Po(a+ •••a+a^i •••a^J ^ Sr,sdet[po{a+aq/)]^^^^^^^ . (3.18) 

That is, any higher order correlation function decomposes into the determinant 
of the matrix of pair correlations. In its most general form, a particle number 
conserving quasifree state po : ^ —>■ C can be written as 

poiA) := ^Tr(e-^A) (3.19) 



for A G a, with 



and 



Zk ■■= Tr(e-«), (3.20) 

K = I dpdqn{p,q)ap Qq (3.21) 

bilinear in a^^aq] for a proof, see [3j. We assume K to be deterministic (with 
respect to {ijJx\x)- 

If in addition, translation invariancc is imposed, such that 

[K,T] = (3.22) 

then 



K ^ dph{p)a+ap (3.23) 

is bilinear and diagonal in n't ,ap. 

Since H^ is bilinear in the creation- and annihilation operators, it is immediately 
clear that 

K{t) := e'*"" Ke~''"" (3.24) 

is also bilinear in a^.aq. Therefore, 

Pt{A) = ^Tr(e-^(*)A) (3.25) 
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is quasifree with probability 1. However, since quasifreeness is a nonlinear condition 
on determinants, almost sure quasifreeness does not imply that E[pt{ ■ )] is quasifree. 

In fact, E[/9j( • )] is not quasifree for any t > 0. 



However, we prove in Theorem l3. 21 below that it possesses a kinetic scaling limit 
(in the sense of Theorem 13. ip which is quasifree. 

Theorem 3.2. Assume that po is number conserving and quasifree, and translation 
invariant. Then, the following holds. For any normal ordered monomial in creation- 
and annihilation operators, 

a+{fi)---a+{fr)a{gi)---aigr), (3.26) 

with r,s G N and Schwartz class test functions fj,gg G SiT'^), and any T > 0, the 
macroscopic 2r-point function 

ni^'-\f,,...Jr;gu...,gr) (3.27) 

:= lini lim E[pT/rj2{a+{fi) ■ ■ ■ a~^{fr) a{gi) ■ ■ ■ a{gr))] 

77— >0 L— ►oo 

exists and is quasifree, 

nP{f^,...,fr;9i,...,9r) =det[4'^(/,,5,)]i<,,^.<,. (3.28) 

The macroscopic 2-point function is the same as in Theorem \3.1\ 

^?{f;9)^fdpFT{p)J{p)9{p), (3.29) 



and Ft{p) solves the linear Boltzmann equation Ii3.12\) with initial condition iS.lSp . 

We note that the assumption of translation invariance can easily be dropped. 
However, we do not address inhomogenous Fermi gases in this text. 



4. Proof of Theorem 13.11 
The proof of Thcorcm l3.1l is obtained from an extension of the analysis in [7lll4j. 

4.1. Duhamel expansion. We consider the Heisenberg evolution of the creation- 
and annihilation operators. We define 

ap{t) := e''"- ape-''"- , (4.1) 

and 

a(/,i) := e'*«-a(/)e-'*«-. (4.2) 

We make the key observation that 

a{f,t) = aift) (4.3) 

where /( is the solution of the 1-particle random Schrodinger equation 

idtft = H^'^ ft := Aft + TjVS'^ft (4.4) 

with initial condition 

/o-/. (4.5) 
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Here, A denotes the nearest neighbor Laplacian on A^, and ifi — H^\^^ is the 1- 
particle Anderson Hamiltonian at weak disorders studied in [7l[8l[T4]. Vu, = K;|ffi 
is the 1-particle multipUcation operator (Kj f){x) = ujxf{x). 

To prove (|4.4p . (|4.5p . we observe that since H^^ is bilinear in a+, a, it follows that 
a{f,t) is a linear superposition of annihilation operators. Therefore, there exists a 
function ft such that a[f,i) — a{ft). In particular, 

idtaUt) = [i^^,a(/t)] 



dp Mp) E{p) ap + -q dp du ft{p)ui{u ~ p) a^ 

= a(A/t) + a(r?l/ji)/,), (4.6) 

and moreover, it is clear that a(/, 0) — a(/o) ~ «(/)• This implies (|4.4|) . (14. 5p . 

Thus, 

Pt(«+(/)a(5)) - po(a+(/,)a(5*)) 

dp dq po ( ftp a<j ) /t (p) 9t{q) 



^ dpJ{p)ft[p)gt{p) (4.7) 

where 

Po{a+aq) = 5{p-q)J{p) (4.8) 

due to translation invariance, with 

< J{p) = ^Poia^ap) = ^^1^ < 1, (4.9) 

cf. p.7p : see p. 231) for the definition of h{p). In particular, this implies (|3.8p . 

For A^ G N, which we determine later, we expand ft, gt into the truncated 
Duhamel series at level A^, 

ft = /i^"^ + /f"^ , (4.10) 

with 

N 



/P'^'-E/* ' (4-11) 

and where the Duhamel term of n-th order (in powers of r/) is given by 

/n 
dso---ds„S{t-Y.'i) (4-12) 

3=0 

/n n 

dko--- dkn s{p - fco) ( n e*^^^("^^ ) ( n '^(% - %-i) ) /(^-) 

== r/"e'* f dae'*" f dko ■ ■ ■ dk^Sip - ko) (4.13) 



■ft _j tL 

n^(fc;y^^7^)(n-(^.-V0)/(M. 

J=0 ^ ■'■' 3 = 1 
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The remainder term is given by 



We choose 



/(>^)=*W d.e^(*-)^^y«/W(.). (4.14) 

Jo 



e = \ (4.15) 



so that the factor e*^* remains bounded for all t. Accordingly, 

N+l 

pt{a+{f)aig)) = po{ a+ (ft) a{gt) ) = ^ pi"'"^(/;s) (4.16) 

n,n—0 

where 

P^'H/;5) :=Po(a+(/("V(gf))) (4.17) 

if n, n < N, and 

pi"'^^^^(/;5) - Po(a+(/i"^)a(5p^^)), 
Pr^'"^(/;5) - Po(a+(/f^V(<7f^)) (4.18) 

if n < A^, respectively if n < A, and 

pr'-^'^'Hf^g) ■■- Po(a+(/P"V(5p"^)). (4.19) 

In particular, for n,n < A, 



dfco • • • *n / d/eo • • • d/cs /(fcn) 9{kn) J(fco) (5(fco - /co) ) 

n -J n -J 

11 "p7r7i~r77~r77 11 



,0 ^(fcj) - a - Je ^^^^p E{ki) - 5 + 



ze 



]^ S(% - %_i) J]^ (:D(fcf_i - h) . (4.20) 

j=i fci 

This expression, and the expressions involving n and / or n = A + 1 , are completely 
analogous to those appearing in the truncated Duhamel expansion of the Wigner 
transform in [71 [M] . 

This permits us to use the methods of [71 [13] to prove Theorem l3.1l We will here 
only sketch the strategy; for the detailed proof, we refer to [71[T3]. In our subsequent 
discussion, we will compare the expressions appearing in the given problem to those 
treated in [3 [14]. 

To begin with, we introduce a more convenient notation. Clearly, if n,n < N, 
and n + n is odd, E[/3( '"'(p, q)] ~ 0. Thus, we let 

n:=^eN, (4.21) 
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and we define {mj}^"o^ by 

kn-j a j < n 
kj-n-i if j > n+ 1 . 

Consequently, 

E[p^"H/;5)] - r;2«e2^* / da d5e**("-") (4.23) 



(4.22) 



duo--- du2n+l /(mo) 9{u2n+l) J{Un) S{Un - U„+l) 



, 2n+l 



Y\ i n 

j=o ^ ■" e=n+i ^ ^' 

n 2fi+l 

e[[|S(wj-Uj_i) [] S(wj-Uj_i)] (4.24) 

j=l j=n+2 

in these new variables, where we use that u){u)* = uj(—u). 

4.2. Graph expansion. Next, we take the expectation with respect to the random 
potential. To this end, we introduce the set of Feynman graphs T„^n, with n + n € 
2N, as follows. 

We consider two horizontal solid lines, which we refer to as particle lines, joined 
by a distinguished vertex which we refer to as the po-vertex (corresponding to 
the term pa{a'^ '^u„+i )• On the line on its left, we introduce n vertices, and on 
the line on its right, we insert n vertices. We refer to those vertices as interaction 
vertices, and enumerate them from 1 to 2n starting from the left. The edges between 
the interaction vertices are referred to as propagator lines. We label them by the 
momentum variables uq, ..., U2n+i, increasingly indexed starting from the left. To 
the j-th propagator line, we associate the resolvent ^r^.-}_^_^^ if < j < n, and 
-gr, — } ~ , - if n + 1 < 7 < 2n + 1. To the £-th interaction vertex (adjacent to the 

edges labeled by U£_i and ui), we associate the random potential Lo{ut, — U£_i), 
where 1 < £ <2n + I. 

A contraction graph associated to the above pair of particle lines joined by the 
Po-vertex, and decorated by n + ri interaction vertices, is the graph obtained by 
pairwise connecting interaction vertices by dashed contraction lines. We denote 
the set of all such contraction graphs by r„ j^; it contains 

(2nV 
|r„,s| = (2fi - l)(2n - 3) • • • 3 • 1 = ^ = Oin\) (4.25) 

elements. 

If in a given graph tt e r„ ;j, the ^-th and the i'-th vertex are joined by a 
contraction line, we write 

i r^^ i' , (4.26) 

and we associate the delta distribution 

6{ui -ue-i - (ue' - U£'-i)) = E[uj{ue - U£-i)Q{u£' - ui'-i)] (4.27) 
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to this contraction line. 



4.3. Classification of graphs. For the proof of Thcorcm l3.H we classify Feynman 
graphs as follows; see [Hill], and Figure 1. 

• A subgraph consisting of one propagator line adjacent to a pair of vertices 
(, and ^ + 1, and a contraction line connecting them, i.e., (. ^t^ ^ + 1, where 
both £, £ + I are either < n or > n + 1, is called an immediate recollision. 

• The graph tt e r„.„ (i.e., n = n ~ n) with £ ~^ 2n — £ for all £ — 1, . . . , ti, 
is called a basic ladder diagram. The contraction lines are called rungs of 
the ladder. We note that a rung contraction always has the form £ ~^ £' 
with £ < n and £' > n+ I. Moreover, in a basic ladder diagram one always 
has that if £i ~^ £[ and £2 ~^ £'2 with €1 < £2, then ^ < £[. 

• A diagram it G r„^fj is called a decorated ladder if any contraction is either 
an immediate recollision, or a rung contraction £j ~^ £', with £j < n and 
£j > n for j = 1, . . . , /c, and £1 <■■■< £k, £[>■■■ > £'k- Evidently, a basic 
ladder diagram is the special case of a decorated ladder which contains no 
immediate recollisions (so that necessarily, n — n). 

• A diagram tt G r„ j^ is called crossing if there is a pair of contractions 
(■ ^T^ ^', J ~ir j', with ^ < f and j < j', such that £ < j. 

• A diagram n G r„^s is called nesting if there is a subdiagram with ^ '^t^ 
£ + 2k, with fc > 1, and either £ > n + 1 or £ + 2k < n, with j ~7r J + 1 for 
j — £+l,£ + 3,...,£ + 2k — 1. The latter corresponds to a progression of 
fc — 1 immediate recollisions. 

We note that any diagram that is not a decorated ladder contains at least a crossing 
or a nesting subdiagram. 



n+1 

1> 



po 



Figure 1. An example of a Feynman graph, n G r„^s, with n = 4, ?i = 6. The 
distinguished vertex is the po-vertex. 



4.4. Feynman amplitudes. Next, we average (|4.20p with respect to the random 
potential. Accordingly, E[ J^tD(M£ — ui-i) ] splits into the sum of all possible prod- 
ucts of pair correlations, according to Wick's theorem (we recall that {ujx} are 
assumed to be i.i.d. Gaussian). This implies that 

E[p|"^")(/;5)] = ^ Amp^(/;5;e;^) (4.28) 

Trer„ s 



2n+1 



with 
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Amp^(/;g;e;77) := yy^n ^2^* f da da e'*^"'-°'^ (4.29) 

duo-- ■ dU2n+l /(wq) giu2n+l) J [u^) 5{Un - U„+i) 
n - 2n ^ 



-'- i?(M,) — a — ie -'■-'- E{ui)—a + ie 



and e = 7- Here, 



-^xlK}-:^^') := W 5{ui-ut-i-{ut,-ut-i)) (4.30) 

is the product of the delta distributions associated to aU contraction hues in tt. 
Moreover, we recall that 

5{ u„ - -u„+i ) J(u„) = po( a+„a«„+i ) , (4-31) 

see (|T5)) . We note that 

UO - U2n+1 = 0, (4.32) 

as one easily sees by summing up the arguments of all delta distributions. This 
holds for any n, n and again implies 



We observe that the role of (|4.3ip in (|4.29p is analogous to that of the rescaled 
Schwartz class function Jg in [ZlIIl], and that the test functions /, g here correspond 
to the initial state (pQ in [7l [14] . 

4.5. Contribution from crossing and nesting diagrams. The amplitude of 
any graph tt G r„ s that contains either a crossing or a nesting can be estimated 

by 

lim |Amp,(/;g;e;,7)| < \\ I h\\ ghW J WooC^" {\og-)\c^\-Hog-r , (4.33) 
L— ►oo e e 

see [71 [13] • We note that similarly as in [71 [T3] , the bounds on all error terms will 
only depend on the L^-norm of the initial condition, which in [7] I14j is normalized 

by ll^olll = 1. 

The existence of the the thermodynamic limit, as L -^ oo, is obtained precisely 
in the same manner as in [T] [S] . Let 

K7 C r„,jj (4.34) 

denote the subset of diagrams of crossing or nesting type. The number of graphs 
in 

r2r := U Kl (4.35) 

n+n— 2n 

is bounded by 2^n!. 
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Thus, the sum of ampHtudes associated to all crossing and nesting diagrams is 
bounded by 

J2 E r^™ |Amp,(/;g;e;77)| (4.36) 



< {N + l)\e^/^\og-)\cr,\-Hog-y 



e e' 

noting that evidently, ||/||2, II3II2 < C for f,g of Schwartz class, and recalling from 
(l3Jl) that 



II J Hoc < 1, (4.37) 

which in particular is the case for J{p) = (1 + e'^^^^P^"''')"^ associated to a Gibbs 
state of the free Fermi field, for all < /3 < 00. 

4.6. Remainder term and time partitioning. If at least one of the indices n, n 

equals A^ + 1 , we first use 

|E[p(^+i-)(/;.g)]| < (E[pf-)(.g;.g)])V2(E[p(^+i^^+i)(/;/)])V2 (4.38) 

by the Schwarz inequality (assuming without any loss of generality that n = A^+1). 
If ri < A^, the term E[p^'^'{g;g)] admits a bound of the form (|4.48p below. 

To bound K[pl ' {f',f)], corresponding to the remainder term in the 
Duhamcl expansion, we use the time partitioning method of [14) : see also [?]• To 
this end, we further expand the remainder term into 3A'^ additional Duhamel terms, 
and to subdivide the time integration interval [0, t] into k G N equal segments 

(4.39) 



(4.40) 
(4.41) 

(4.42) 
&^''^^\P) ■■= i^vT-'' / dso--- dsn-N S{ V s, - {s - t)) (4.43) 

/n-N n-N 

duQ--- dUn-N S{p - Uq) JI e'^J-EK) "Q q(^^_ _ ^^._^^ /(u„-7v) ■ 
i=0 £=1 

Moreover, 







wherel 


yy one obtains 




A>N) ^ ANAN) ^ .(>4JV) 


where 


K AN-l 
j=l ,i=Ar+l 


with 


Jin,N,r,..) :=.r7yW/i-^--^), 


and 
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We note that writing (I4.4ip in the form J2'j=i J2n=N+i 5" J' ^^ have 

K 4N-1 

- J2 Yl ko(a^(5n',i')a(ff«j))| 

_7 j' — 1 n,n' —N-\-l 
K 4Ar-l 

j^j' — l n.7i'—N-\-l 

< k'^ (SN)"^ suppo(a^(gnj)a(g„j)). 
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(4.44) 



Thus, by the Schwarz incquahty, 



^(^+i,Ar+i)^^.^^ <i?l(/,t) + i?2(/,t) 



(4.45) 



where 



i?i(/,t) := (3iV)2«2 sup po(a+(4"'^'--^))a(/("^^^--^))) (4.46) 



Af<Ti<47V 
1<3<K 



and 

i?2(/,t):=i^ sup sup po(a+(/1^^4^'--^))a(/1^^4^---^))). (4.47) 

l<i<Kse[rj_i,rj] 

By separating terms due to decorated ladders from those due to crossing and nesting 
diagrams, one finds 



hmE[po(a+(/4"'^'--^))a(/(;-^---^)))] 

= E[y'dpj(p)|/^;^^^--^)(p)|2; 



(4.48) 



< \\i\\i\\J\ 



^f0J^+.^/^(log^)^(c,V4ogi)- 



for A^ < n < AN (see [71 [H] for a detailed discussion) . 



For n — AN, the main issue is to control the large factor i^ in (|4.47p . To this 
end, we observe that for a time integral on the interval [tj^i,tj] of length -, the 
parameter e = t~^ can be replaced by ne = (-)~^. Therefore, one gets 



limE[po(a+(/r'"^''^-^^)a(/r'"^'^^-^^))] 



< ||J||oo]E[||/]^'^^'^-^Ml 

< II/II2IIJII 



«^^)^logi)^c.^.-logl)- 



.2N 



(4.49) 



The gain of a factor k ^^ is crucial; it is sufficient to compensate for the factor P 
in (|4.47p , using the parameter choice given in Section 14.71 below. 
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One obtains that if at least one of the indices n, n equals iV + 1, 

£m^|E[/"-) (/;/)] I < H/||^HJ||oo[ '' (jV!)i/2 

+ (tV^ «:2 gl/5 + ,-2 ^-2Ar ^ ((4^),) (log i)4(^;^2^-l ^^g i )8Ar j ^ (4 ^q) 

where k remains to be chosen. The first term on the right hand side of (|4.50p 
bounds the contribution from all basic ladder diagrams contained in the Duhamel 
expanded remainder term. For a detailed discussion, we refer to [71I51I14J. 



4.7. Choosing the constants. We recall from ()4.37|) that || J ||oo < 1- Moreover, 
II / II2, II g II2 < C for aU test functions /, g e S{T'^). As in [71[1[I1], we choose 



1 T 



e rf 



N 



10 log log i 



(log-)i^ (4.51) 



e 



Then, 



and consequently, 



and 



e-1/11 < ATI < e-1/10 

K^ > £-3/2 (4.52) 

(14361) , (144911 < 77^/1^ (4.53) 



(113011 < 77I/4 (4.54) 



for 77 sufficiently small. It follows that the sum of all crossing, nesting, and remain- 
der terms is bounded by 77^/^°. 



4.8. Resummation of decorated ladder diagrams. Let T^^"'-' C r„ ^ denote 
the subset of all decorated ladders based on n + n vertices. Then, for T > 0, let 

W(£(T,,,)) 

^P''\f;9) ■■= J2 J2 J2 ^lim Amp,(/;g;e(T,77);77) (455) 

with e{T,ri) = \- In the kinetic scaling limit 77 ^ with t = \ = T/rf, one 
obtains 

4')(/;.g) :- lim ^'^''^/j.g) - [ dp Ft{p)J{p) g{p) , (456) 

7,^0 J 
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where 



Ft{p) := \iinF!f\p) 



27tT J du{E{u)-E{p)) 



(4.57) 



J2 f^^^^^dSo---dSnS{T~J2^^) 

n=0 "'K+^^ j=0 

/ft 
duo ■ ■■ dunS{p-uo)( ]^2 7r(5(-E(wj) -£;(wj_i)) J Fo(u„) 



with initial condition 



Fo(u) = Km J{uA-) = lim 7^ Po(a„^. au^* ) 



(4.58) 



(for the definition of ua* , see Theorem 13. ip . It can be straightforwardly verified 
that (|4.57p is a solution of the Cauchy problem for the linear Boltzmann equation 
(|3.12p . as asserted in Theorem 13. II 



5. Proof of Theorem EH 



Because both K (in the definition of po) f^nd the random Hamiltonian H^^ are 
bilinear in a"^, a (of the form J dui du2 k{ui, U2) a^ a^j), the same is true for 



K{t) := e "^ K e 



with probability 1. Therefore, 

Pt(-) = ^Tr(e-^(*)(.)) 
is quasifree with probability 1 (see, for instance, [3j). Thus, for r, s e N, 



/Ot(a^(/i)---a^(/r)a(5i)---a(.9s)) ^Sr^sdet pt{a^ {Jj)a{gt)) 



3,i=l 



(5.1) 



(5.2) 



(5.3) 



where fj,gi G iS(T'^) belong to the Schwartz class. In particular, we can set r = s. 
We expand the determinant into 



det 



Pt{a'^{f])a{gt)) 



j-A=i 



(5.4) 



where Sr is the symmetric group of degree r. We claim that for T > and t — ^^ 
and any choice of fj^gt G iS(T'^), 



lim 

L — >oo 



^[\{pT/n^{a+{fj)a{g,(j)))'\^ - J| E[pT/r,2(a+(/j) a(5,(,)) ] | < r;^(5.5) 
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for a constant S > independent of r, s e 5r, t], and T, and for 77 > sufficiently 
small. This immediately implies that, for every fixed r < cxd, 



lim 

L — ^oc 



E 



PT/,,2{a+{fi) ■ • •a+(/r) a(.gi) • • • a{gr) ) 



det 



¥.[pT/r,2ia+{fj)a{g(,) 



j.fci 



< r! 77 



converges to zero as 77 ^ 0. 
This implies that 

J2r 



^T ih,---Jr;gi,---,gr) 



lini lim E[pt/ 2 ( a+(/i) • ■ • a+(/,.) 0(51) ■ • • 0(5,.) )] 



is quasifree, i.e. 

l(2r) 



^T (A'---'/'-; gi,---,9r) = det[ny{fi; gj)]^ 



^i,j^r 



where 



>(2)r 



f^F(/;5) = dpFT{p)f{p)g{p). 



(5.6) 



(5.7) 



(5.8) 



(5.9) 



The function Ft{p) solves the linear Boltzmann equation with initial condition 
-Fo(p), as given in Theorem 13. II 



5.1. Proof of (|5.5p . The inequality (|5.5p follows from a straightforward applica- 
tion of the main results in [5] where we refer for details. In this section, we shall 
only outline the strategy. The expectation 



lim ¥.\\{pt{a+{f,)a{g,^,))) 

j=i 



(5.10) 



can be represented by a graph expansion as follows. We expand each of the factors 

(5.11) 



Pt{a+{fj)a{gs(j))) = 



N+l 

E 

n,n— 1 



f(") 



(H) 



dpj{p)ff;\p)g';^l^^^{p) 



separately into a truncated Duhamel series of level N ^ using the same definitions 
as in (|4.16|) . For the remainder term (where at least one of the indices n,n equals 
N + 1), we subdivide the time integration interval [0, t] into k, pieces of length -. 



For the expectation (|5.10p . we introduce the following extension of the classes of 
Feynman graphs discussed for the proof of Theorem l3.11 see also Figure 2. For r > 1, 
we consider r particle lines parallel to one another, each containing a distinguished 
Po-vertex separating it into a left and a right part. Enumerating them from 1 to 
r, the j-th particle line contains rij interaction vertices on the left of the po-vertex, 
and rij interaction vertices on its right. We note that for r > 1, only X]j=i('^i +^i) 
has to be an even number, but not each individual 

rij := rij + rij . (5-12) 

On the j-th interaction line, we label the propagator lines by momentum variables 



,U) 



,(i) 



i- , , , with indices increasing from the 



left. 
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^ O 



n^+1 



po 



n-i+1 



^ c> 



^ ^ 



n2+1 



po 



^j o 



n2+1 



^ ^ 



"3 ■ "3+^ ^ 



a — I 1> 



po 



n^+l 



^ ^ 



n +1 

r 



o o 



po 



n+1 



Figure 2. Order r Feynman graph. The particle line indexed by j = 3 is discon- 
nected. 



A contraction graph of degree {{nj,nj)}^^-^ is obtained by connecting pairs of 
interaction vertices by contraction lines. We denote the set of contraction graphs 
of degree {(fij,?^j)}j=i by T^(^nj,nj)}':_ ■ If the £-th vertex on the j-th particle line 
is connected by a contraction line to the ^'-th vertex on the j'-th particle line, we 
write 



U-J) -. {f ;(■')■ 

To a graph tt G '^{(n-,n ■)}':_ 1 we associate the Feynman amplitude 

Amp,({/„g,(,)};r;;r) := r;2 2:,<,<.K+",) ^2... f[ J da, da, e^'^^^-"- 



(5.13) 



,0) ...^„0') f.{.,(Jh 



,{^'> \.l(-„i3)\S(nM) -„(^) 



d<' ■ --dy^Ui /.(<0 .9sO)«/+i) Jiul^pS{ul^> uli>^, ) (5.14) 

1 t\ 1 



sAb^^}Z^')I[ 



n 



E{u'f^) - aj - ie ^,^„.+2 E{u'f?) - a, + ie 



where 



1 rj^ 
'=1 = T 



(5.15) 
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for T > 0. The delta distribution 

5AWpr;^)= n 5{v!f^-v^\-iuf:''-v^:i^) (5.16) 

(3;i) ~. U'-i') 

is the product of deha distributions associated to all contraction lines in tt. 



5.1.1. Completely disconnected graphs. The subclass 

^fin,,n,)r,^, C r{(„^.^jj,)}^-^, (5.17) 

of completely disconnected graphs of degree {('T-j,%)}r=i consists of those graphs 
in which contraction lines only connect interaction vertices on the same particle 
line. 

It is clear that 



lii^ E E Amp,({/„5,(,)};r;;r) (5.18) 

r N „ 



i— 1 nj ,nj — 1 



= lim n (E[pT/,^(a+(/,)a(g,(,)))] + 0(77^)) , (5.19) 

i=i 

according to our proof of Theorem 13.11 The term of order 0{irj^) accounts for 
the remainder term associated to the j'-th particle line (i.e., the terms involving 
E[p^y"" (p, g)] where at least one of the indices nj,nj equals N). Thus, for any 
fixed r S N, we obtain 

lim lim V V Amp^({/j,.g^(j)};77;r) 

Tj >{j 1j — * CXID 

j = l,...,r -K"-j>"j)>j = i 

r 

= n4"(/.;5s0-)). (5.20) 

That is, the sum over completely disconnected graphs yields the corresponding 
product of averaged 2-point functions in the kinetic scaling limit. 

5.1.2. Non-disconnected graphs. We refer to the complement of the set of com- 
pletely disconnected graphs in T^i^^ ,n )}^_ , 

^{K-,S,)}-^, — r|(„^,K,)}J^, \ ^{(n],n,)}-^,^ (5.21) 

as the set of nan- disconnected graphs. It remains to prove that the sum over non- 
disconnected graphs, combined with the remainder terms, can be bounded by 0(77 ), 
for L sufficiently large. 

The condition required in [8] for the estimate analogous to (|5.5p to hold is that 
for the initial condition 4>q (corresponding to the test functions fj , gi in our case) 
of the random Schrodinger evolution studied in [8], a "concentration of singularity 
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condition" is satisfied (that is, singularities in momentum space are not too much 
"spread out" in the limit 77 — > 0). It states that in frequency space T"*, 

00 = $t'^ + $1:'"'^ , (5.22) 

where 

II ^('■^^) lU < c (5.23) 

and 

|||-(..„,)|^|-(..„,)|||^ <c'^3/^ (5.24) 

are satisfied uniformly in L, as L ^ cx). 

In the present case, we have to require that fj, gi satisfy the concentration 
of singularity condition. This is, however, evidently fulfilled since fj, gi are rj- 
independent Schwartz class functions (in contrast, the initial states considered in 
[8] are of WKB type, and scale non-trivially with rj.) 

It is proven in 8J that the amplitude of every non-disconnected graph with 
rij, Hj < N for J = 1, . . . , r, is bounded by 

sup I Amp^({/j,5r3(j)};77;r) I (5.25) 

< ei/5(c772e~i logi)5 2:j=i».(iogi)4^, (5.26) 

e e 

where we recall that ^ = 7 = 7^ for T > 0. This key estimate is a factor e^/^ 
smaller than the bound on the sum of disconnected graphs; this improvement is 
obtained from exploiting the existence of at least one contraction line that connects 
two different particle lines; see [8]. 

The number of non-disconnected graphs is bounded by 

|r^(-4,.,)}^J < (E%)! < (2-^)! (5-27) 

where Uj = rij + rij. Therefore, the sum of amplitudes of all non-disconnected 
graphs with < n^ , rij < N is bounded by 

E E E |Amp,({/„.g,(,)};r;;T)| (5.28) 

< ((2r7V)!)2ei/5 (^^2 ^-i ^^g If ^ {log i)^'- . (5.29) 

Here we have estimated the sum over pairs < rij, rij < iV, 1 < j < r, by another 



factor (2riV)!, since #{inj,nj)}'j^^ \ J2j^j = "i} < 



m 



5.1.3. Duhamel remainder term. In case at least one of the indices rij or rij equals 
N + 1, the following argument can be applied. Clearly, from a Holder estimate of 
the form \\hi ■ ■ ■ hs\\i < \\hi\\s ■ ■ ■ \\hs\\s with respect to E, we have 

E[f[pi-^'^^\f;g)]\ < nE[|pi"-^^)(/;.g)p'-]-- (5.30) 
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Here, we have used an exponent 2r instead of r because then, even for r odd, an 
absolute value of the form jzp'' can be replaced by a product of the form 'z^z^ , 
where z £ C. 

We make a choice of constants 

1 T 

t = - = 



e rf 



N 



log 7 



lOrloglogi 

(logi)i5-, (5.31) 

e 



similarly as in Section 14.71 of the proof of Theorem 13.11 

If rij or rij equals iV + 1, we can use the bounds (|4.53[) and (I4.54p . 
If both n, n < N, we use the a priori bound 

Y^ ^ limE[|p("'")(/;g)pn* (5.32) 



_ _ L—*OG 

n+n— 2n TV^Tn^n 

2r 



^2^~l\n.i 



<[E(T)(^W7^| ^ 



e e 



X e^/^ ((2r - e)ny. ( (log -f{cTjh-^ log i)" V 



< ^|i]T7f + '^^ (5-33) 

The factor ,^^y — in [■•■] accounts for i basic ladders on £ copies of ^^Ji^ 
while the remaining factor accounts for all other (not necessarily non-disconnected) 
contractions on the remaining 2r — £ particle lines; for details, see [71I51I14J. 

Let us without any loss of generality assume that ni = A^ + 1. Then, keeping ni 
fixed and summing over the remaining indices ni and rij, rij, with j = 2, . . . , r, we 
find 

r 



y: n^[ipr'"^-^(/;5) 



0<n2,nj,nj<N + l j = l 
j = 2,...,r 



<^?"iE \Ui/2 +^"J (5-34) 

n=0 ^ '' 

where the factor 77 ts accounts for the remainder term indexed by rii = A'^ + 1. We 
conclude that the sum over all terms (|5.30p which contain at least one rij or rij 
equalling iV + 1 (i.e., which contain at least one Duhamcl remainder term) can be 
bounded by 

Crjrs (5.35) 

for a constant C independent of 77 and r. 
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Combined with 

({5:28)) < rj^ , (5.36) 



which one easily verifies, this completes the proof of Theorem l3.2l For more details 
addressing the arguments outlined here, we refer to [5j . 
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